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ABSTRACT 



A search for diffuse neutrinos with energies in excess of 10^ GeV is conducted 
with AMANDA-II data recorded between 2000 and 2002. Above lO'^ GeV, the 
Earth is essentially opaque to neutrinos. This fact, combined with the limited 
overburden of the AMANDA-II detector (roughly 1.5 km), concentrates these 
ultra high-energy neutrinos at the horizon. The primary background for this 
analysis is bundles of downgoing, high-energy muons from the interaction of cos- 
mic rays in the atmosphere. No statistically significant excess above the expected 
background is seen in the data, and an upper limit is set on the diffuse all-flavor 
neutrino flux of $9o%cl < 2.7 x 10^'' GeV cm^^ s^^ sr^^ valid over the energy 
range of 2 x 10^ GeV to 10^ GeV. A number of models which predict neutrino 
fluxes from active galactic nuclei are excluded at the 90% confidence level. 

Subject headings: neutrino telescope, AMANDA, IceCube, diffuse sources, ultra 
high-energy 



1. Introduction 



AMANDA-II (Antarctic Muon and Neutrino Detector Array), a neutrino telescope 
at the geographical South Pole designed to detect Cherenkov light from secondary parti- 
cles produced in collisions between neutrinos and Antarctic ice, h as placed limits on the 
flux from point-like and diffuse sources of astrophysical neutrinos (lAchterberg et al.l 120071 : 
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Ackermann et al.ll2004j . l2005l : lAhrens et al.ll2003al Jbl). This work describes a search for neu- 
trinos with energies above 10^ GeV, which we define as uhra high-energy (UHE) neutri- 
nos. These neutrinos are of interest because the y could be associated with the potential 



acceleration of hadrons by active galactic nuclei (IMannheinj Il995l : iMannheim et al.l 12000 



Halzen fc Zas 1997 : Protheroe 1996 : Stecker et al. 19921 ): they cou ld potentially be produced 
by exotic phenomena such as the decay of topological defects (ISigl et al.l Il998l ) or possi- 
bly associated with the Z-burst mechanism (lYoshida et al.lll998l ): and they are guaranteed 
by-prod ucts of the interac tions of high-energy cosmic rays with the cosmic microwave back- 
ground (lEngel et al.ll200ll ). 



This analysis is sensitive to all three fiavors of neutrinos. Leptons and cascades from 
UHE electron, muon and tau neutrinos create bright, energetic events (Fig. [T]) which can 
be identified by AMANDA-II as far as 450 m from the center of the array (Fig. [2]). The 
sensitivity of this analysis starts at energies roughly coincident with the highest energy 
threshold of other diff use analyses conducted with AMANDA-II (lAchterberg et al.l 120071 : 
Ackermann et al.ll2004l ). 



At UHE energies, the interac tion length of neutrinos in rock is shorter than the diameter 
of the Earth (iGandhi et al.lll998l ). so neutrinos from the Northern Hemisphere will interact 
before reaching AMANDA-II. Combined with the limited overburden above the AMANDA- 
II detector, this concentrates UHE events at the horizon. This contrasts with the majority of 
other astrophysical neutrino analyses completed using data from the AMANDA-II detector, 
which search for neutrinos from the Northern Hemisphere with energies below 10^ GeV. 

The fiux of atmospheric neutrinos is negligi ble at UHE e nergies, with fewer than 10 
events in three years expected from the model in iLiparil (119931 ) after intermediate UHE se- 
lection criteria have been applied. This drops to 0.1 events after application of all selection 
criteria. Similarly, there are fewer than 0.6 events expected in three years at the final se- 
lection level from prompt neutrinos fron i the decay of ch armed particles produced in the 
atmosphere (using the "C" model from IZas et al.l (119931 )). Therefore, the primary back- 
ground for the UHE analysis is composed of many lower energy processes that mimic higher 
energy signal events. Cosmic ray collisions in the upper atmosphere that generate large 
numbers of nearly parallel muons (or "muon bundles") can generate high-energy signatures 
even though the individual muons have much lower energy than single leptons or cascades 
from UHE neutrinos. Signal and background events spread light over roughly equivalent 
areas in the detector, but UHE neutrino events are distinguishable because they have higher 
energy and higher light density than background events. Specialized selection criteria which 
use these properties, as well as differences in reconstruction variables, separate the UHE 
neutrinos from the background of muon bundles from atmospheric cosmic rays. 
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Limits have been placed on the all-flavor neutrino flux in the ultra high-energy range 
by other experiments (Fig. [3]). Additionally, a previous analysis using an earlier configu- 
ration of the AMANDA detector called AMANDA-BIO, consisting of 302 optical modules 
( lAckermann et al.ll2005l ). has placed limits on the all-flavor UHE neutrino flux (Fig. [3]). This 
analysis uses 677 optical modules (OMs) of the AMANDA-II detector and gives a combined 
result using data from three years (2000-2002) with a livetime of 456.8 days. 

A description of the AMANDA-II detector is given in section 2. Sections 3 and 4 discuss 
possible sources of astrophysical neutrinos and background, and the simulation of both. The 
selection criteria used to separate UHE neutrino signals from background are discussed in 
section 5. A study of systematic uncertainties is presented in section 6, and the results are 
shown in section 7. 



The AMANDA-II Detector 



The AMANDA-II detector ( lAhrens et al.ll2004al ) consists of 677 OMs stationed between 
1500 m and 2000 m beneath the surface of the Antarctic ice at the geographic South Pole. 
The OMs are deployed on nineteen vertical cables (called "strings" ) arranged in three roughly 
concentric circles, giving the detector a cylindrical shape with a diameter of approximately 
200 m. 

Each OM contains a Hamamatsu 8-inch photomultiplier tube (PMT) coupled with sil- 
icon gel to a spherical glass pressure housing for continuity of the index of refraction. The 
OMs are connected to the surface by cables which supply high voltage and carry the sig- 
nal from the PMT to data acquisition electronics at the surface. The inner ten strings use 
electrical anal og signal transmissio n, while the outer nine strings primarily use optical fiber 
transmission (lAhrens et al.ll2004al ). 



The AMANDA-II detector uses a majority trigger of 24 OMs recording a voltage above 
a set threshold (a "hit") within a time window of 2.5 fis. An OM records the maximum 
amplitude, as well as the leading edge time and time over threshold for each hit, with each 
OM recording a maximum of eight hits per event. Each photoelectron has approximately a 
3% chance of prod ucing an afterpulse caused by ionization of residual gas inside the PMT 
(IHamamatsul Il999l ) . This afterpulse follows several /is after the generating hit and aids in 
the detection of UHE events. 

AMANDA-II has been collecting data since February 2000. In 2002/ 2003, waveform 
digitizers were installed which record t he full pulse shape f rom each OM (jSilvestril 120051 ). 
In 2005 deployment began on IceCube (lAhrens et al.ll2004bl ). a 1 km^ array of digital OMs 
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which now encompasses the AMANDA-II detector. 



3. Astrophysical Neutrino and Background Sources 



Astrophysical neutrinos with energies in excess of 10^ GeV may be produced by a 
variety of sources. A number of theories predict neutrino fluxes from active galactic nuclei 
(AGN) peaking near 10^ GeV. In these scenarios, protons are accelerated by the first order 
Fermi mechanism in shock fronts. In the favored mechanism for neutrino production , these 



Protheroe 


1996; 


Halzen & Zas 


1997: 



+ fj^ ^ + + Ue + v^,, 



and produce neutrinos via the process: 

p + 7 ^ A+ ^ 7r+[+n] 

resulting in a Uf.-.u^-.i'T- fiavor ratio of 1:2:0 at the sourc^ll- The energy spectrum of the 
neutrinos produced by these interactions generally follows the spectrum of the protons. 
Theoretical bounds can be placed on the flux of these neutrinos based on the observation 
of cos mic rays if the p-7 reaction takes place in the jet or other optically thin region of the 
AGN jBahcall fc WaxmanlEggsl : iMannheim et allboooh . 



UHE neutrinos are also associated with models created to explain the apparent excess 
of cosmic rays at the highest energies. One scenario involves the decay of massi ve objects, 
such as topological defects created by symmetry breaking in the early universe (jSigl et al. 
19981 ). These objects decay close to the Earth into showers of particles, eventually producing 
neutrinos as well as a fraction of the highest-energy cosmic rays. Z-burst models could 
also produce some of the highest-energy cosmic rays through the interaction of neutrinos 
with energies in excess of 10^^ GeV with relic neutrinos via the TP resonance. Since these 
neutrino- neutrino interactions are rare, it is possible to directly search for the UHE neutrino 
fluxes required by this mechanism (lYoshida et al.lll998l : iKalashev et al.ll2002al ). It should be 
noted that Z-burst scenarios which pr edict the highest flu x of neutrinos have already been 
eliminated by previous experiments (IBarwick et al.l l2006l ). Additionally, Z-burst models 
predict fluxes of neutrinos which peak at energies above the sensitivity of this analysis or 
require unrealistic assumptions and are mentioned primarily for completeness. 

A guaranteed source of UHE neutrinos comes fro m the intera c tion o f hig h-energy cosmic 



rays with the cosmic microwave background (see e.g. lEngel et al.l (1200 ll ) and lKalashev et al. 



^Neutrino flavor oscillation changes the flavor ratio to 1:1:1 at the Earth. See Kashti fc Waxman ( 2005[ ) 
for a discussion of different flavor ratios. 



- 8 - 



( l2002bl )). However, the flux predictions of tliese GZK neutrinos are generally several orders 



of magnitude lower than most of the fluxes listed previously. 

The background for this analysis consists of bundles of muons from cosmic rays. The 
cosmic rays foll ow an E"^ '' spectrum until about 10^ GeV, where the flux steepens to E^^ 



( lHorandelll2003l ) . They come only from the Southern Hemisphere because bundles from other 
directions are absorbed by the Earth. According to simulations, there can be as many as 
20,000 muons in one bundle spread over a rms cross-sectional large as 200 m^, and 

the highest-energy events can deposit energies as large as 2.4 x 10^ GeV in the ice around 
the AMANDA-II detector. 



4. Simulation and Experimental Data 



UH E neutrinos are simulated using the All Neutrino Interaction Simulation (ANIS) 
package (iKowalski &: Gazizovi l2005l ) to generate and propagate the neutrinos through the 
Earth. All three flavors of neutrinos are simulated with energies between 10^ GeV and 
10^^ GeV. The resulting muons and taus are propagated through the rock and ice near 
the d etector using the Muon Monte Carlo (MMC) simulation package (jChirkin &: Rhode 
20041). Finallv. the d etector response is simulated using the AMASIM2 simulation package 
( iHundertmarkl 1 1 998l ) . 



The background muon bundles from cosmic rays are generated using the CORSIKA 
simulation program with the QGSJETOl hadronic interaction model (jHecklll999l ). At early 
levels of this analysis, cosmic ray prim aries are generated with composition and spectral 
indices from IWiebel-Sooth et al.l (119991 ). with energies of the primary particles ranging be- 
tween 8 X 10^ GeV and 10^^ GeV. At later levels of this analysis, the lower energy primaries 
have been removed by the selection criteria, and a new simulated data set is used with en- 
ergy, spectral shape, and composition optimized to simulate high-energy cosmic rays more 
efficiently. In this optimized simulation, the energy threshold is raised to 8 x 10^ GeV and 
only proton and iron primaries are generate d with a spectrun a of E~ ^. These primaries are 
reweighted following the method outlined in iGlasstetter et al.l (119991 ). This optimized simu- 
lation is used for level 2 of the analysis and beyond (see Tabled]). For 2001 and 2002, the 
background simulation is further supplemented with the inclusion of a third set of simulated 
data with the energy threshold increased to 10^ GeV. For all sets of background simulation, 
the resulting particles are propagated through the ice using MMC, and the detector response 
is simulated using AMASIM2. 



Data used in this analysis were recorded in the time period between February 2000 and 
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November 2002, with breaks each austral summer for detector maintenance, engineering, 
and cahbration lasting approximately four months. In addition to maintenance downtime, 
the detector also has a brief period while recording each event in which it cannot record 
new events. Runs with anomalous characteristics (such as excessive trigger rates or large 
numbers of OMs not functioning) are discarded and a method which remo ves non-physical 
events caused by short term detector instabilities is applied (jPohJ 1200J) . These factors 
combine to give a deadtime of 17% of the data taking time for 2000, 22% of the data-taking 
time for 2001, and 15% of the data taking time for 2002. Additionally, 26 days are excluded 
from 2000 because the UHE filtered events are polluted with high number of events with 
incomplete hit information, likely due to a minor detector malfunction. Taking these factors 
into account, there are 173.5 days of livetime in 2000, 192.5 days of livetime in 2001, and 
205.0 days of livetime in 2002. Finally, 20% of the data from each year is set aside for 
comparison with simulations and to aid in the determination of selection criteria, leading to 
a total livetime for the three years of 456.8 days. 



5. Analysis 



Twenty percent of the data from 2000 to 2002 (randomly selected from throughout the 
three years) is used to test the agreement between background simulations and observations. 
In order to avoid biasing the determination of selection criteria, this 20% is then discarded, 
and the developed selection criteria are applied to the remaining 80% of the data. A previous 
UHE analysis was perforr ned on ori l y the 200 data using diff erent selection criteria than 
those described below (see lGerhardtl (120051 ) and lGerhardtl (120071 ) for a more detailed descrip- 
tion). For 2001-2002 , improved reconstruction techniques such as cascade reconstructions 
(lAhrens et al.ll2003d ) were added to the analysis, and the new selection criteria described 
below were devised in a blind manner. These selection criteria were also applied to the 2000 
data to derive a combined three year limit. Due to differences in hit selection for reconstruc- 
tion between 2000 and 2001-2002, the signal passing rate at the final selection level for 
the year 2000 is approximately 60% of the rate for the years 2001 and 2002. 

In order to maximize the limit setting potenti al, the selection crite ria are initially de- 
termined by optimizing the model rejection factor (IHill fc RawlinsI 120031 ) given by 



MRF 



N, 



signal 



(2) 



where figo is 90% confidence level (CL) average event upper limit given by lFeldman &: Cousins 
( 119981 ) ■ and Ngignai is the number of muon neutrinos expected for the signal being tested, in 
this case an E~^ flux. The selection criteria for this analysis are summarized in Table [1] and 
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described below. 

This analysis exploits the differences in total energy and light deposition between bun- 
dles of many low-energy muons and single UHE muons or cascades from UHE neutrinos. 
UHE neutrinos deposit equal or greater amounts of light in the ice than background muon 
bundles. In addition to being lower energy, background muon bundles spread their light 
over the cross sectional area of the entire muon bundle, rather than just along a single muon 
track or into a single cascade. Both signal and background events can have a large number 
of hits in the array, but for the same number of hit OMs, the muon bundle has a lower total 
number of hits, NHITS (recall each OM may have multiple separate hits in one event; see 
Fig. H]). The number of hits for UHE neutrinos is increased by the tendency of bright signals 
to produce afterpulses in the PMT. Background muon bundles also have a higher fraction of 
OMs with a single hit (FIH), while a UHE neutrino generates more multiple hits (Fig. 
The FIH variable is correlated with energy (Fig. [6]) and is effective at removing lower energy 
background muon bundle events. The level 1 and 2 selection criteria require that NHITS > 
140 and FIH < 0.53 and reduce the background by a factor of 2 x 10'^ relative to trigger 
level (level on Tables [2] and [3]) . 

At this point the data sample is sufficiently reduced that computationally intensive 
reconstructions become feasible. Reconstruction algorithms used in this analysis employ 
a maximum likelihood method which takes into account the absorption and scattering of 
light in ice. For muons, the reconstruction cor apares time residuals to those expected from a 



Cherenkov cone for a minimally ionizing muon (lAhrens et al.ll2004al ) , while the c ascade recon- 



struct ion uses Cherenkov light from an electromagnetic cascade for comparison (lAhrens et al. 



2003d ). Reconstructions which are optimized for spherical (cascade) depositions of light are 
used to distinguish UHE neutrinos from background muon bundles which happen to have a 
large energy deposition, such as a bremsstrahlung or e"*" e~ pair creation, inside the detector 
fiducial volume. 

Before application of the level 3 selection criteria, the data sets are split into "cascade- 
like" and "muon-like" subsets. This selection is performed using the negative log likelihood 
of the cascade reconstruction (Lease, see Fig. [7]), where events with a Lease < 7 are considered 
"cascade-like." 



5.1. "Cascade-like" Events 



Background events in the "cascade-like" subset are characterized by either a large light 
deposition in or very near the instrumented volume of AMANDA-II or a path which clips 



-li- 



the top or bottom of the array. In either case, the energy deposition is significantly less 
than the energy deposited by a UHE neutrino, allowing application of selection criteria 
which correlate with energy. One of these is FIHelec (Fig- ED, a variable similar to the 
FIH variable described above, except that it uses only OMs whose signal is brought to the 
surface by electrical cables. The signal spreads as it propagates up the cable, causing hits 
close together in time to be combined. This gives FIHelec a different distribution from 
FIH, and both variables are good estimators of energy deposited inside the detector (Fig. 
E]). Additionally, the fraction of OMs with exactly four hits (F4H) is another useful energy 
indicator. The value of four hits was chosen as a compromise between the number of hits 
expected from OMs with electrical cables and OMs with optical fibers. OMs with optical 
fibers typically have more hits than OMs with electrical cables because very little pulse 
spreading occurs as the signal propagates up the fiber. The level 3 selection criteria uses the 
output of a neural net with FIHelec, F4H, and FIH as input variables (Fig. [H]). As selection 
levels 4 and 5, separate applications of the F4H and FIHelec variables remove persistent 
lower energy background events. 

The remaining background muon bundles have a different hit distribution than UHE 
neutrinos. In the background muon bundles, a large hght deposition can be washed out 
by the continuous, dimmer light deposition from hundreds to tens of thousands of muons 
tracks. In contrast, UHE muons can have one light deposition that is several orders of 
magnitude brighter than the light from the rest of the muon track and looks very similar 
to bright cascades from UHE electron and tau neutrinos. For all cases, the initial cascade 
reconstruction is generally concentric with this large energy deposition, so ignoring OMs 
that are within 60 m of the initial cascade reconstruction reduces the fraction of OMs that 
are triggered with photons from the cascade. For background, the remaining light will be 
dominated by light depositions from the tracks of the muon bundles and be less likely to 
reconstruct as a cascade. In contrast, signal events, with their energetic cascades, will still 
appear cascade- like and reconstruct with a better likelihood (Lgo)- The final selection criteria 
for "cascade-like" events (chosen by optimizing the MRF) requires that these events be well 
reconstructed by a cascade reconstruction performed using only OMs with distances greater 
than 60 m and reduces the background expectation to events for this subset. 

The number of events at each selection level for experiment, background, and signal 
simulation for the "cascade-like" subset are shown in Table [2l 
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5.2. "Muon-like" Events 

Background events in the "muon-like" subset are characterized by more uniform, track- 
hke hght deposition and are more easily reconstructed by existing reconstruction algorithms 
than "cascade-like" events. A reconstruction algorithm based on parameterization of time 
residuals from simulated muon bundles is used to reconstruct the zenith angle of the events 
(Fig. [TOj) . Since most background muon bundles will come from a dow ngoing direction, whil e 



UHE neutrinos will come primarily from the horizontal direction (IKlein fc MannI Il999l ). 
requiring that the zenith angle > 85° (where a zenith angle of 90° is horizontal) reduces 
the background by a factor of 30. The remaining background in the "muon-like" subset 
are misreconstructed events, since the actual flux close to the horizon is very small. A 
reconstruction based on the hit pattern of a Che renkov cone for a minimally ionizing muon 



is applied to these events (lAhrens et al.l l2004a] ). Selecting only well-reconstructed events 



using the hkelihood of this reconstruction (Lmuon) is sufficient to remove all background 
events in this subset. The value of this selection criteria was initially chosen to optimize 
the MRF for muon neutrinos with an E^^ spectrum. However, by increasing the selection 
value slightly beyond the value which gave the minimum MRF, all background events were 
rejected with only a few percent drop in the sensitivity (Fig. [TTI) . Since the uncertainty in 
the cosmic ray spectrum is very large at these energies, the more stringent selection criterion 
was applied to correct for the fact that the MRF is optimized without uncertainties. 

The number of events at each selection level for experiment, background, and signal 
simulation for the "muon-like" subset are shown in Table [31 



6. Statistical and Systematic Uncertainties 



Because there is no test beam which can be used to determine the absolute sensitivity 
of the AMANDA-II detector, calculations of sensitivity rely on simulation. The dominant 
sources of statistical and systematic uncertainty in this calculation are described below. 
The systematic uncertainties are assumed to have a fiat distribution and are summed in 
quadrature separately for background and sig nal. The uncertainties have been included into 
the final limit using the method described in iTegenfeldt fc Conrad J2005h . 



6.1. Uncertainties Due to Limited Simulation Statistics 



Due to computational requirements, background simulation statistics are somewhat lim- 
ited. Ideally, one would scale the statistical uncertainty on zero events based on the simula- 
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tion event weights in nearby non-zero bins. However, the optimized background simulations 
used in this analysis have large variations in event weights approaching this region, making 
determination of this factor difficult. Nevertheless, the statistical uncertainties near the edge 
of the distribution are on the order of the uncertainties for a simulation with a livetime equiv- 
alent to the data taking period, so no scaling factor is applied to the statistical uncertainty. 
A stati stical uncertainty of 1.29, t he la Feldman-Cousins event upper limit on zero observed 
events (IFeldman fc Cousins! Il998l ). is assumed at the final selection level. Signal simulation 
has an average statistical uncertainty of 5% for each neutrino fiavor. 



6.2. Normalization of Cosmic Ray Flux 



The average energy of cosmic ray primaries at the penultimate selection level is 4.4 
X 10^ GeV, which is considerably above the knee in the all-particle cosmic ray spectrum. 
Numerous experiments have measured a large spread i n the absolute normalization of the 
fiux of cosmic rays at this energy (see iKampertI (120071 ) for a recent review ). Estimates o f 
the uncertainty in t he normalization of the cosm ic ray fiux range from 20% (jHorandelll2003l ) 
to a factor of two (jParticle Data Groupl l2006l ). This analysis uses the more conservative 
uncertainty of a factor of two. 



6.3. Cosmic Ray Composition 



There is considerable uncertainty in the cosmic ray composition above the knee (jParticle Data Group 



20061 ) . We estimate the systematic uncertainty by considering two cases: proton-dominated 
composition and iron- dominated composition. The simulated background cosmic ray fiux 
is approximated by separately treating proton and iron primaries combined in a total spec- 
tr um that becomes effecti vely iron-dominated above 10^ GeV using the method described 



m 



Glasstetter et al.l (119991 ). The iron-dominated spectrum yields a 30% higher background 



event rate than the rate from a proton-dominated spectrum at the penultimate selection 
level. This value of 30% is used as the uncertainty due to the cosmic ray composition. 



6.4. Detector Sensitivity 



The properties of the refrozen ice around each OM, the absolute sensitivity of individual 
OMs, and obscuration of OMs by nearb y power cables can e ffect the detector sensitivity. 
This analysis uses the value obtained in lAhrens et al.l (j2003al ) where reasonable variations 
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of these parameters in the simulation were found to cause a 15% variation in the E ^ signal 
and background passing rate. 



6.5. Implementation of Ice Properties 



As photons travel through the ice they are scattered and absorbed. The absorption 
and scattering lengths of the ice arou nd the AMANDA-II det ector have been measured very 
accurately using in situ light sources (lAckermann et al.ll2006l ). Uncertainties are introduced 
due to the limited precision with which these parameters are included in the simulation. 
Varying the scattering and absorption lengths in the detector simulation by 10% were found 
to cause a difference in number of expected signal events (for an spectrum) of 34% 



(lAckermann et al.ll2005l ). which is used as a conservative estimate of the uncertainty due 
to implementation of ice properties. If too large of a deviation in background rate relative 
to the experimental rate was observed for a set of ice property parameters, the background 
rate was normalized to the experimental rate, and the signal rate was scaled accordingly. 
This was done to ensure that the variation in absorption and scattering lengths covered a 
reasonable range of ice properties. 



6.6. Neutrino Cross Section 



T he uncertainty in the stan dard model neutrino cross section has been quantified re- 
cently (lAnchordoqui et al.ll2006l ) taking into accou nt the experimerital un certainties on the 
parton distribution functions measured at HERA (jChekanov et al.ll2005l ). as well as theo- 
retical uncertainties in the effect of heavy quark masses on the parton distribution function 
evolution and on the calculation of the structure functions. The corresponding maximum 
variation in the number o f expected signal events (for an E~^ spectrum) is 10%, in agreement 
with previous estimates (lAckermann et al.ll2005l ). 



Screening effects are expec ted to suppress the n e utrino - nucleon cross s ection at energies 
in excess of 10* GeV (see e.g. iKutak fc Kwiecihskil (120031 ): iBerger et al.l (120071 )). This has 
a negligible effect on the number of signal events expected for an E~^ spectrum because 
the majority of signal is found below these energies f Fig. [T2|). Even if the suppressi on is 



as extreme as in the Colour Glass Condensate model (iHenley &: Jalilian-Marian 
event rate decreases by only 11%. 



2005f ). the 
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6.7. Differences in Simulated Distributions 

An examination of tlie Lmuon distribution for tlie "muon-like" subset after level 3 of this 
analysis suggests the background simulation is shifted by one bin relative to the experiment 
(Fig. [T3l) . Shifting all simulation distributions to the left by one bin leads to better agreement 
between the background simulation and experimental distributions and an increase in 8% in 
the number of expected signal events for an spectrum. 



6.8. The Landau-Pomeranchuk-Migdal (LPM) Effect 



At ultra high-energies, the LPM effect suppresses the bremsstrahlung cross section for 
electrons and the pair-production cross section of phot ons created in a cascade by an electron 
neutrino (ILandau &: Pomeranchuklll953l : iMigdall 119571 ). This lengthens the resultant shower 
produced by a factor that goes as ^/E. Above 10^ Ge V, the exten ded shower length becomes 
comparable to the spacing between OMs on a string (jKleinll2004l ). Additionally, as the LPM 
effect suppresses the bremsstrahlung and pair productions cross sections, photonuclear and 
electronuclear interactions begin to dominate which lead to the production of muons inside 
the electromagnetic cascade. Toy simulations were performed which superimposed a muon 
with an energy of 10^ GeV onto a cascade with energy of 10^ GeV. While the addition of 
the muon shifted the Lease distribution 5% towards higher (more "muon-like") values, the 
resulting events still passed all selection criteria indicating that the effects of muons created 
inside cascades are negligible. 

The LPM effect is not included in the simulations of electron neutrinos, but it can be 
approximated by excluding all electron neutrinos with energies in excess of 10^ GeV. This 
is an overestimation of the uncertainty introduced by the LPM effect, as extended showers 
may manifest as several separate showers which are likely to survive all selection criteria 
and the addition of low-energy muons is not expected to significantly alter the UHE cascade 
light deposition. Neglecting electron neutrinos with energies in excess of 10^ GeV reduces 
the number of expected signal events by 2% for an spectrum. 



6.9. Summary of Uncertainties 

The systematic errors are shown in Table IH Summing the systematic errors of the signal 
simulation in quadrature gives a systematic uncertainty of ±39%. Combining this with the 
statistical uncertainty of 5% per neutrino flavor gives a total maximum uncertainty of 40%. 
Following a similar method for the background simulation, the systematic uncertainty is 
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+ 101% / -60%. Scaling the statistical uncertainty of the background simulation by the 
systematic uncertainty gives a maximum background expectation of fewer than 2.6 events 
for three years. 



Results 



After applying all selection criteria, no background events are expected for 456.8 days. 
Incorporating the statistical and systematic uncertainties, the background is expected to 
be found with a uniform prior probability between and 2.6 events. A possible sensitivity 
calculation which incorporates these uncertainties can be generated by assuming a flat prior 
with a mean of 1.3 events and a corresponding data expec tation of 1 event. This gives a 
90% CL event upper limit of 3.5 (ITegenfeldt fc Conradll2005l ) and a sensitivity of 1.8 x 10~^ 
GeV cm~^ s~^ sr~^, with the central 90% of the signal found in the energy range 2 x 
10^ GeV to 10^ GeV. Table shows the expected number of each flavor of UHE neutrino 
passing the final selection level for a 10^^ x E^^ flux. The energy spectra of each flavor are 
shown in Fig. 



Two events are observed in the data sample at the flnal selection level (Fig. [T3l) . while 
fewer than 2.6 background events are expected which gives a 90% CL average event upper 
limit of 5.3. After applying all selection criteria, 20 events are expected for a 10~^ x E^^ 
all flavor flux (Table [5]). The upper limit on the all-flavor neutrino flux (assuming a 1:1:1 



: Ut flavor ratio) is 



E^'^'90%CL < 2.7 X lO^^GeV cm^^ 



sr 



(3) 



including systematic uncertainties, with the central 90% of the E ^ signal found between the 
energies of 2 x 10^ GeV and 10^ GeV. 

A number of theories which predict fluxes with non-E~^ spectral shapes (Fig. [3]) were 
al so tested by reweight ing the simulated signal. These include the hidden-core AGN model 



of IStecker et al.l ( 



emission (IStecker 



99 



^ which has been updated to reflect a better understanding of AGN 



2005h. as well as AGN models in whi c h neutrinos are accelerated in op 



ticall y thin regions (jProtherod Il996l : iHalzen fc Zad 119971 : iMannheinj Il995l : iMannheim et al. 



2000h. Includ i ng un cert ainties, this analysis res tricts at a 90% CL the AGN models from 



20051 ) models from 



Protherod ( 



Halzen fc Zad (11997) andlMannheim et al.l (120001). Also the prev iously rejected (lAckermann et al 



1996h and Istecker et all (Il992|) are rejected at the 90% CL 



by this analysis (see Fig. UM and Table [6]). The model by IStecker et al.l (119921 ) builds on a 
correlation between X-rays and neutrinos from AGNs. Other models using the same corre- 
lation give a similar normalization and violate current limits by an order of magnitude as 
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well. As previously pointed out by lBecker et al.l (120071 ). such a correlation can be excluded 



While we do not directly exclude the flux from the ISteckerl (120051 ) hidden-core AGN 
model, it is possible to set limits on the parameters used in the model. In this model, the 
flux of neutrinos is normalized to the extragalactic MeV photon flux measured by COMPTEL 
with the assumption that the flux of photons from Seyfert galaxies is responsible for 10% 
of this MeV background. If the neutrino flux scales linearly with the photon flux, then the 
maximum contribution of hidden-core AGNs, such as Seyfert galaxies, to the extragalactic 
MeV photon flux must be less than 29%. 



Fluxes of neutrinos from the decay of t opological defects (|Sigl et al.lll998r) and the UHE 
fluxes required for the Z-bursts mechanism ( Yoshida et al. 1998r Kalashev et al. 2002a) peak 
at too high of an energy to be detected by this analysis. Neutrinos from the interaction of 
cosmic rays with cosmic microwave background photons are produced at too low of a flux 
for this analysis to detect (see Table [6]). 

The number of expected events of a given flavor (z/ and V) for spectra not tested in this 
paper can be calculated using the formula 



signal 



(4) 



where T is the total livetime (456.8 days), A^^ is the angle averaged neutrino effective area 
(Fig. [T5l) . and is the flux at the Earth's surface. 



8. Conclusion 



The diffuse neutrino flux limit for a 1:1:1 : u^j^ : v^- flavor ratio set by this analysis of 

E^*90%cL < 2.7 X lO^^GeV cm'^ s"^ sr~\ (5) 

is the most stringent to date above 10^ GeV. A number of models for neutrino production 
have been rejected (see Table |6] for a full list). AMANDA-II hardware upgrades which were 
completed in 2 003 should lead to an improvement of the sensitivity at ultra high-energies 
( ISilvestril |2005| ) . Additionally, AMANDA-II is now surrounded by the next-generation Ice- 
Cube detector which is currently under construction. The sensitivity to UHE muon neutrinos 
for 1 year is expected to increase by roughly an order o f magnitude as the IceCube detector 
approaches its final size of 1 km^ (lAhrens et al.ll2004bl ). 



We acknowledge the support from the following agencies: National Science Foundation- 
Office of Polar Program; National Science Foundation-Physics Division; University of Wis- 



- 18 - 



consin Alumni Research Foundation; Department of Energy, and National Energy Research 
Scientific Computing Center (supported by the Office of Energy Research of the Depart- 
ment of Energy); the NSF-supported TeraGrid system at the San Diego Supercomputer 
Center (SDSC); the National Center for Supercomputing Applications (NCSA); Swedish 
Research Council; Swedish Polar Research Secretariat; Knut and Alice Wallenberg Founda- 
tion, Sweden; German Ministry for Education and Research; Deutsche Forschungsgemein- 
schaft (DFG), Germany; Fund for Scientific Research (FNRS-FWO); Flanders Institute to 
encourage scientific and technological research in industry (IWT); Belgian Federal Office for 
Scientific, Technical and Cultural affairs (OSTC); the Netherlands Organisation for Scien- 
tific Research (NWO); M. Ribordy acknowledges the support of the SNF (Switzerland); A. 
Kappes and J. D. Zornoza acknowledges the Marie Curie OIF Program; L. Gerhardt ac- 
knowledges the support of the University of California, Irvine MPC Computational Cluster 
and Achievement Rewards for College Scientists (ARCS). 



REFERENCES 

Achterberg, A. et al., 2007 Phys. Rev. D, in press flarXiv:0705.1315p . 

Ackermann, M. et al., 2004, Astroparticle Physics, 22, 127. 

Ackermann, M. et al., 2005, Astroparticle Physics, 22, 339. 

Ackermann, M. et al., 2006, Journal of Geophysical Research, 3, D13203. 

Ahlers, M. et al., 2005, Phys. Rev. D, 72, 023001. 

Ahrens, J. et al., 2003, Astrophys. Journal, 583, 1040. 

Ahrens, J. et al., 2003, Phys. Rev. Lett., 90, 251101. 

Ahrens, J. et al., 2003, Phys. Rev. D, 67, 012003. 

Ahrens, J. et al., 2004, Nuclear Instruments and Methods A, 524, 169. 

Ahrens, J. et al., 2004, Astroparticle Physics, 20, 507. 

Anchordoqui, L., Cooper-Sakar, A., Hooper D., and, Sakar, S., 2006, Phys, Rev. D, 74, 
043008. 

Aynutdinov, V. et al., 2006, Astroparticle Physics, 140. 
Bahcall, J. and Waxman, E., 1998, Phys. Rev. D, 59, 023002. 



- 19 - 



Barwick, S. et al., 2006, Phys. Rev. Lett., 96, 171101. 

Becker, J. K., Grofi, A., Munich, K., Dreyer, J., Rhode, W. and Biermann, P. L., 2007, 
Astroparticle Physics, in press ( astro-ph/0607427 ) . 

Berger, E., Block, M., McKay, D., and Tan, C, 2007, preprint (hep-ph/0708.1960vl). 

Chekanov, S. et al., 2005, European Physical Journal C, 42, 1. 



Chirkin, D. and Rhode, W., 2004, preprint ( |hep-ph/0407075| ). 

Engel, R., Seckel, D., and Stanev T., 2001, Phys. Rev. D, 64, 093010. Curve shown for 
A=0.7, taken froni |ft^:/7ftp.bartol.udel.edu/seckel/ess-gzk/flux_n3_8Jlat_oni0p3.txt| 

Feldman, G. and Cousins, F., 1998, Phys. Rev. D, 57, 3873. 

Gandhi, R., Quigg, C, Reno, M., and Sarcevic, I., 1998, Phys. Rev. D, 58, 093009. 



Gerhardt, L., 2005, Proc. 29th Int. Cosmic Ray Conf., Pune, India, 111 ( astro-ph/0509330) . 

Gerhardt, L., 2007, Proc. 14th Int. Conf. on Supersymmetry and the Unification of Funda- 
mental Interactions, AIP Conf. Proc. 903, 622. 

Glasstetter, R. et al., 1999, Proc. 26th Int. Cosmic Ray Conf., Utah, USA, 1, 222. 

Halzen, F. and Zas, E., 1997, Astrophys. Journal, 488, 669. 

Hamamatsu, 1999, Photomultiplier Tubes, Basics and Applications, second ed. 

Heck, D., 1999, DESY-PROC-1999-01, 227. 

Henley, E. and Jalilian-Marian, J., 2005, preprint ( |hep-ph/0512220| v^l). 
Hill, G. and Rawlins, K., 2003, Astroparticle Physics, 19, 393. 
Horandel, J., 2003, Astroparticle Physics, 19, 193. 

Hundertmark, S., 1998, Proc. 1*** Workshop Methodical Aspects of Underwater/Ice Neutrino 
Telescopes, Zeuthen, Germany. 

Kalashev, O. et al, 2002, Phys. Rev. D 65, 103003. 

Kalashev, O. et al., 2002, Phys. Rev. D 66, 063004. 

Kampert, K.-H., 2007, Nuclear Physics B (Proc. Suppl.), 165, 294. 



- 20 - 



Kashti, T. and Waxman, E. Phys. Rev. Lett. 95, 181101. 
Klein, J. and Mann, A., 1999, Astroparticle Physics 10, 321. 
Klein, S., 2004, preprint ( |astro-ph/04l2546l vl). 

Kowalski, M. and Gazizov, A., 2005, Computer Physics Communications, 172, 203. 
Kravchenko, I. et al., 2006, Phys. Rev. D, 73, 082002. 



Kutak, K. and Kwiecihski, J., 2003, preprint ( hep-ph/0 303209y4). 
Landau, L. and. Pomeranchuk, I., 1953, Dok. Akad. Nauk SSSR, 92, 535. 
Lipari, P., 1993, Astroparticle Physics, 1, 195. 
Mannheim, K., 1995, Astroparticle Physics, 3, 295. 

Mannheim, K., Protheroe, R. J., and Rachen, J., 2000, Phys. Rev. D, 63, 023003. 
Migdal, A., 1957, JETP, 5, 527. 

Particle Data Group, 2006, Journal of Phys. G. 33, 1. 
Pohl A.C., 2004, Licenciate thesis, Kalmar University. 



Protheroe, R., 1996, preprint ( astro- ph/96071 65). 

Silvestri, A., 2005, Proc. 29th Int. Cosmic Ray Conf., Pune, India, 431. 
Stecker, F.,2005, Phys. Rev. D, 72, 107301. 

Stecker, F., Done, C, Salamon, M., and Sommers, P., 1992, Phys. Rev. Lett., 69, 2738. 

Sigl, G., Lee, S., Bhattacharjee, P., and Yoshida, S., 1998, Phys. Rev. D, 59, 043504. 

Tegenfeldt, F. and Conrad, J., 2005, Nuclear Instruments and Methods in Physics Research 
A, 539, 407. 

Wiebel-Sooth, B., Biermann, P., and Meyer, H., 1999, col. VI/3c, Springer Verlag, 37. 

Yoshida, S., Sigl, G., and Lee S., 1998, Phys. Rev. Lett., 81, 5505. 

Zas, E., Halzen, F., and Vazquez, R., 1993, Astroparticle Physics, 1, 297. 



This preprint was prepared with the AAS I^TJrjX macros v5.2. 



-21 - 



Table 1: Selection criteria. 



Level 


Selection Criteria 





Hit Cleaning and Retriggering 


1 


FIH < 0.72 




NHITS > 140 


2 


FIH < 0.53 


"Cascade-like" "Muon-like" 


3 


7^7 T > 7 

-'-'case ^ ' -'-'case _ ' 




Neural Net > 0.93 Zenith Angle > 85 


4 


F4H < 0.1 Wuon < 6.9 


5 


FIHelec < 0.56 


6 


Leo < 6.6 
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Table 2: Number of experimental, simulated background, and simulated signal events in the 
"cascade-like" subset at each selection level for 456.8 days. 



Level Data BG Simulation Signal Simulation 






2.7 xlO^ 


1.8^};^ xlO^ 


621.7 


1 


3.9 xlO^ 


3.1+?] xlO^ 


270.8 


2 


1.7 xlO^ 


1-4^^1 xlO^ 


89.2 


3 


155 


62+p 


32.0 


4 


151 


f5i+62 
Di_37 


31.0 


5 


46 


qr)+32 


27.1 


6 





0+2.6 


16.0 



Note. — Levels and 1 show combined numbers for both "muon-hke" and "cascade-hke" subsets. 
Signal is shown with a low energy threshold of 10'* GeV for a neutrino spectrum of dN/dE = 10~^ x 
E^^ GeV^^ cm^^ s^^ sr^*, with an assumed 1:1:1 : i/^ : flavor ratio. Values at selection level and 
1 for data and background simulation are extrapolated from the 2000 datasets. The background simulation 
is shown with systematic and statistical uncertainties described in Section [6l The number of "muon-like" 
events are shown in Tabled 
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Table 3: Number of experimental, simulated background, and simulated signal events in the 
"muon-like" subset at each selection level for 456.8 days. 



Level Data BG Simulation Signal Simulation 






2.7 xlO^ 


i.8tl:l xio' 


621.7 


1 


3.9 xlO^ 




270.8 


2 


1.4 xlO^ 


Q.Ol^i xlO^ 


85.2 


3 


4.6 xlO^ 


2.7tl-l xlO^ 


57.9 


4 


2 


0+2.6 


4.0 



Note. — Levels and 1 show combined numbers for both "muon-hke" and "cascade-hke" subsets. 
Signal is shown with a low energy threshold of 10^ GeV for a neutrino spectrum of dN/dE = 10~^ x 
E^^ GcV^^ cm^^ s^^ sr~^, with an assumed 1:1:1 Ve v^j. Vt flavor ratio. Values at selection level and 
1 for data and background simulation are extrapolated from the 2000 datasets. The background simulation 
is shown with systematic and statistical uncertainties described in Section [6l The number of "cascade-like" 
events are shown in Tabled 
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Table 4: Simulation Uncertainties 



Source 



BG Sim 



Sig Sim 



Cosmic Ray Normalization +100% / -50% 

Cosmic Ray Composition -30% 

Detector Sensitivity ±15% ±15% 

Ice Properties - ±34% 

Neutrino Cross Section - ±10% 

Simulation Distribution - ±8% 

LPM Effect - -2% 



Total 



+101% / -60% +39% / -39% 
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Table 5: Number of simulated neutrino events in the "cascade-like" and "muon-like" subsets 
passing all selection criteria for three years for a neutrino spectrum of d(N,^^-|-N,^^-|-N,^^)/dE 
= 10-6 X E-2 GeV -1 cm-2 s'^ sr'^ 



Neutrino Flavor 


"Cascade-like" 


"Muon-hke" 


Total 


Electron 


7.7 


0.1 


7.8 


Muon 


3.9 


3.6 


7.5 


Tau 


4.4 


0.3 


4.7 


All Flavors 


16.0 


4.0 


20.0 
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Table 6: Flux models, the number of neutrinos of all flavors expected at the Earth at the 
final selection level, and the MRFs for 456.8 days of livetime. 



Model 




MRF 


Reference 


AGN^ 


20.6 


0.3 


fProtheroe 1996) 


AGN'^ 


17.4 


0.3 


fStecker et al. 1992) 


AGN'^ 


8.8 


0.6 


fHalzen & Zas 1997) 


AGN'^ 


5.9 


0.9 


f Mannheim et al. 2000) 


AGN RL 


4.5 


1.2 


fMannheim 1995) 


Z-burst 


2.0 


2.7 


fKalashev et al. 2002a) 


AGN 


1.8 


2.9 


fStecker 2005) 


GZK V norm AGASA^ 


1.8 


2.9 


fAhlers et al. 2005) 


GZK V mono-energetic 


1.2 


4.4 


fKalashev et al. 2002b) 


GZK V a=2 


1.1 


4.8 


fKalashev et al. 2002b) 


GZK V norm HiRes*^ 


1.0 


5.3 


fAhlers et al. 2005) 


TD 


0.9 


5.9 


fSid et al. 1998) 


AGN RL A^ 


0.3 


18.0 


fMannheim 1995) 


Z-burst 


0.1 


53.0 


fYoshida et al. 1998) 


GZK V 


0.06 


88.0 


fEnerel et al. 2001) 



"These values have been divided by two to account for neutrino oscillation from a source with an initial 1:2:0 
v^:v^j_:v-r flux. 

^Lower energy threshold of 10'' GeV applied. 

Note. — A MRF of less than one indicates that the model is excluded with 90% confidence. 



e 




Fig. 1. — Simulated muon neutrino event with an energy of 3.8 x 10^ GeV. The muon passes 
roughly 70 m outside the instrumented volume of the detector. Colored circles represent hit 
OMs. The color of the circle indicates the hit time (red is earliest), with multiple colors 
indicating multiple hits in that OM. The size of the circle is correlated with the number of 
photons produced. 
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Muon distance of closest approach [m] 



Fig. 2. — Distance of closest approach to the detector center for muons from UHE muon 
neutrinos (shown with an spectrum) which pass all selection criteria of this analysis. 
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Fig. 3. — All-flavor UHE neutrino flux limit for 2000-2002 over the range which con- 
tains the central 90% of the expected signal with an spectrum. Also shown 
are s everal rep r esenta tive models: St05 from ISteckerl (120051) mul t iplied by 3, P96 
from iProtherod (1l996l ) multiplie d by 3/2, En g Ol fr om lEngel et al. ( 20011 ), Si98 from 



Sigl et al.1 (119981), Yosh98 from lYoshida et all fll998h. LiD93 from Ihiparil (ll993|), and 



the Waxman-Bahcall upper bound f Bahcall &: Waxrnan 19981) multipl i ed by 3/2. Exist- 
i ng experimental lim its shown a re from the RICE (jKravchenko et al.M2006l ). ANITA-lite 
( Barwick et al. 2006), and Baikal (Avnutdinov et al. 20061 ) experiments, the UHE limit from 



A MANDA-BIO (Ackermann et al. 20051 ), the lower-energy diffuse mu on limit multiplied by 



3 (lAchterberg et al.l 120071 ) and cascade limit (lAckermann et al.l 120041 ) from AMANDA-II. 
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Fig. 4. — NHITS distribution for the experiment, background, and E ^ muon neutrino signal 
simulations before level 1 of this analysis. 
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Fig. 5. — Distribution of FIH (the fraction of OMs with a single hit) for the experiment, 
background, and muon neutrino signal simulations after level 1 of this analysis. The 
average FIH drops with energy (see Fig. [6]). 
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Fig. 6. — FIH (top) and FIHelec (bottom) distributions for various energy decades of muon 
neutrino signal. These variables serve as rough estimator of energy for the UHE analysis. 
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Fig. 7. — Distribution of Lease for the experiment, background, and electron, muon, 
and tau neutrino signal simulations after level 2 of this analysis. Events with Lease < 7 are 
"cascade-like," and events with Lease > 7 are "muon-like." 
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Fig. 8. — FIHelec (the fraction of electrical OMs with a single hit) distribution for the 
experiment, background, and electron, muon, and tau neutrino signal simulations in the 
"cascade-like" subset after level two of this analysis. 
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Fig. 9. — Distribution of neural net output for the experiment, background, and electron, 
muon, and tau neutrino signal simulations in the "cascadc-likc" subset after level two of this 
analysis. Signal events are expected near one, while background events are expected near 
zero. 
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Fig. 10. — Reconstructed zenith angle distribution for the experiment, background, and 
electron, muon, and tau neutrino signal simulations in the "muon-like" subset after level two 
of this analysis. Zenith angles of 90° correspond to horizontal events, and zenith angles of 
0° are downgoing events. 
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Fig. 11. — Model rejection factor for 10 ^ x E ^ muon neutrinos in the "muon-like" subset 
as a function of cut level for Lmnon ■ 
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Fig. 12. — Energy spectra of electron, muon, and tau neutrino signal events 
(d(N,,^+N,,^+N,,J/dE = 10~^ x E-^ GeV ~^ cm"^ s"^ sr~^) which pass all selection cri- 
teria. The peak in the electron neutrino spectrum just below 10^ GeV is due to the Glashow 
resonance. 
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Fig. 13. — Lmuon distribution for the experiment, background, and 4.5 x 10~^ x muon 
neutrino signal simulations (arbitrary normalization) in the "muon-like" subset after level 
three of this analysis. Two experimental events survive the final selection criteria of Lmuon 
< 6.9. 
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Fig. 14. — AGN fluxes tested in this paper. Lines denote the model predictions and symbols 
denote the 90% CL hmits on the pred ictions derived by thi s analysis. The model s reject ed at 



the 90% CL s: 
MPROO from 



lown are: HfcZ97 frora iHalzen fc Zad (119971 ). P96 from iProtherod ( 119961 ). and 



Mannheirn et al.l (120001 ). Also shown a r e mo dels close to being rejected: M95 

RL B from Mannheim ( 19951 ) and St05 from Stecker ( 2005 ). See Table |6]for exact numbers. 
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Fig. 15. — Angle- averaged neutrino effective area at final selection level as a function of 
neutrino energy for electron, muon, and tau neutrinos. The peak in the electron neutrino 
effective area just below 10^ GeV is due to the Glashow resonance. 



